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Abstract 

We compare the results of the fixed-flavor scheme calculation of Laenen, Riemersma, Smith and van Neerven 
with the variable-flavor scheme calculation of Aivazis, Collins, Olness and Tung for the case of neutral-current 
(photon-mediated) heavy-flavor (charm and bottom) production. Specifically, we examine the features of both 
calculations throughout phase space and compare the structure function F2(x,Q 2 ). We also analyze the depen- 
dence of F-2. on the mass factorization scale fi. We find that the former is most applicable near threshold, while 
the latter works well for asymptotic Q 2 . The validity of each calculation in the intermediate region is dependent 
upon the x and Q 2 values chosen. 
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Figure 1: The basic semi-inclusive deeply inelastic scattering process, + N(P) — > + Q{pi) + X(Px). 

1 Introduction 

Several experimental groups jlj] have studied the semi-inclusive deeply inelastic scattering (DIS) process for heavy- 
quark production (Figure |l|) 

h{h) + N(P) -> £ 2 (£ 2 ) + Q(pi) + X{P X ). (1) 

Most analyses of this process assume that the hadron is comprised only of the massless gluon (g), up (it), down (d) 
and strange (s) quarks, while the charm (c), bottom (6), and top (t) quarks are treated as massive objects which are 
strictly external to the hadron. 

This view of the heavy quarks as external to the hadron is appropriate when the energy scale of the process /x p h y 
(for example the center-of-mass energy y/s) is not large compared to the mass of the heavy quark, i.e. MQ^$/i p hy. For 
most fixed-target facilities, this condition holds for the c, b and t quarks j|. We are therefore justified in excluding 
c, b, and t as constituents of the hadron in the QCD-improved quark-parton model (QPM) for this case. 

With new data from HERA, the electron-proton collider at DESY, we can investigate the DIS process in a very 
different kinematic range from that available at fixed-target experiments In this new realm, the important 
question is: Should the c and b quarks be considered as partons, or as heavy objects extrinsic to the hadron? Given 
that HERA extends the kinematic reach of the DIS process by two orders of magnitude, we can not expect our 
assumptions that were valid for fixed-target processes to hold in a completely different kinematic regime. 

Aivazis, Collins, Olness and Tung (ACOT) have discussed this issue at length in reference Q and approach the 
problem by invoking the variable flavor scheme (VFS), which varies the number of partons according to the relevant 
energy scale /i p h y - The fundamental physical insight to the VFS is that in the region Mq 3> ^phy, the heavy quark 
should be excluded as a constituent of the hadron as it is kinematically inaccessible. However, when Mq -C Mphy 
the heavy quark should be included as a parton since Mq is insignificant compared to /i p hy Although the physics 
is unambiguous in these kinematic extremes, most experimental data lies in between these clear-cut regions. In 
the intermediate region, the renormalization scheme of Collins, Wilczek and Zee (CWZ) g provides a well-defined 
transition between these two extreme kinematic domains. 

The issue of what constitutes a parton also points to an inconsistency between traditional charged-current and 
neutral-current heavy-quark production calculations Q. When considering charged-current processes, one begins 
with the purely electroweak process W + q — > Q, as shown in Figure @a. For neutral-current processes, the traditional 
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Figure 2: |^a) The O(ot° s ) flavor excitation (FE) process, shown here for the charged-current case, 
flavor creation (FC) process, shown here for the neutral-current case. 



|b) The 0{a\) 
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Figure 3: Pictorial representation of the fixed flavor scheme (FFS) calculation of LRSN. ||a) The 0{a\) flavor 
creation (FC) contribution. ||b), ||c), ||d), and ||e) An illustrative subset of the complete set of 0{a 2 s ) diagrams 



contributing to the complete LRSN calculation. 

approach is to begin with the 0{a\) process 7* + g — > Q + Q, as shown in Figure ||b.[] As we work in the new 
kinematic regime spanned by HERA, the concept of a 'heavy' quark becomes a relative term depending upon the 
magnitudes of the kinematic variables involved. Traditional distinctions between the charged-current and neutral- 
current calculations should vanish as the characteristic energy scale becomes significantly larger than the heavy-quark 
mass. ACOT implements the CWZ renormalization scheme and treats both charged-current and neutral-current 
heavy-flavor production in the consistent fashion of beginning both calculations at 0(a^). 

Laenen, Riemersma, Smith and van Neerven (LRSN) have calculated heavy-quark production for DIS photon 
exchange, beginning with the 0{a\) photon-gluon fusion process and including the complete 0{a 2 s ) radiative correc- 
tions in reference ||. Figure ^| displays some of the relevant 0(g 2 ) tree- level and 0(g 3 ) virtual Feynman diagrams 
involved, where g is the QCD coupling strength and is given by a s — g 2 /Air. LRSN calculates from the viewpoint 
that the produced heavy quark is generated from an initial-state gluon or lighter quark. For example, in producing 
c quarks, LRSN assumes the parton initiating the heavy-quark production process to be the g, u, d or s quark. The 
same tenet holds for b production, only the c quark would be included as a massless initial-state parton as well. 

Another feature of the LRSN calculation is that it provides additional information on inclusive differential cross 
section distributions. At 0(al), the heavy quark and antiquark are produced back to back in the 7*-parton center 
of momentum frame, and at 0(a 2 ), the additional influences of the gluon radiation and the 7*9 channel have also 
been calculated in M. LRSN gives additional insight into the differential distributions, which is particularly useful 
from the experimental point of view. 

^Due to the traditional inconsistency between the neutral and charged processes discussed above, the terminology leading-order and 
next-to- leading-order is ambiguous. Therefore, we denote the subprocesscs according to the power of a s . 
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Figure 4: Pictorial representation of the variable flavor scheme (VFS) calculation of ACOT. |Ja) The O(a s ) flavor 
excitation (FE) contribution. [|b) The schematic representation of the subtraction term. The x on the intermediate 
quark line indicates a convolution between the C( a s) 9 ~ > 1 splitting (lower portion of the diagram) with the 0{oP s ) 
process (upper portion of the diagram). [|c) The 0{a 1 s ) flavor creation (FC) contribution. 

The subject of this paper is a vigorous comparison of the advantages of each calculation, as well as a glimpse 
of future prospects including a merging of the two calculations to produce a three-order result which should have 
excellent predictive power for heavy-flavor structure functions at HERA. 



2 Kinematics 



While Figure [j] shows the general deeply inelastic scattering (DIS) process, we shall focus specifically on neutral- 
current heavy- flavor production via a photon exchange^ as described by the sub-process 7*(<z) + N(P) — ► Q(pq) + 
X(Px)- We define Q 2 = —q 2 , where q is the four momentum of the virtual photon exchanged, y — (P • q)/(P ■ i\) is 
the fractional energy transfer, and x — Q 2 /(2P ■ q) is the Bjorken x variable. The cross section is obtained from the 
structure functions via 

d 2 a 



dx dy 



2MEi G X G 2 
7r m 



9 2 +i 



xFi y + F 2 



*F 3 y(l-l) 



(2) 



Ei is the energy of the incoming lepton, M is the mass of the hadron being probed, Gi = g%./{Q 2 + is a 

shorthand for the boson coupling and the propagator, ng is the number of polarization states of the incoming lepton, 
and <7j_ ; = g 2 Ll ± gj^, where gu and gm refer to the chiral couplings of the vector boson to the leptonsj-*] 



3 ACOT Calculation 

The ACOT calculation makes use of the O(a s ) and 0{a\) processes to obtain a two-order result, cf. Figure || and 
r igure | The O(a s ) result uses the standard QCD evolution to resum the iterative gluon and quark splittings to 
give the parton distribution function (PDF) of the heavy quark. 
The 0(a° s ) term is straightforward, and is represented by 

°F = F(Q(qg)) = /^®-F C0) (7*Qi^ft) , (3) 

* The ACOT calculation with general masses and couplings applies to both charged and neutral-current processes. 

§In this definition of the structure functions, we have extracted the quark coupling (iJ+j) and the average over the incoming lepton 

polarization (ne) so that the same formula applies to both charged and neutral processes. 
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where denotes the parton distribution function for a nucleon N producing a quark Qi, and the F^{^*Qi — > Qi) 
is the lowest-order partonic structure function for the initial-state heavy quark Qi to absorb a virtual photon 7* and 
enter the final state. We introduce the shorthand F(0(a™)) = n F. The <g> represents the convolution which is the 
integral over the momentum fraction carried by the initial state parton i 

fh®F = £ F^.a.^)) ■ (4) 

where £ = k + / P + is the ratio of the "+" -component of the parton momentum fc + to the hadron momentum P + in 
light cone coordinates, /j, r is the renormalization scale, and fj,f is the mass factorization scale. For the remainder of 
this paper, we shall not distinguish between \i r and fif, and we choose [x r = /U/ = fx. Although the partonic structure 
function F and the evolution of the PDF with mass factorization scale \i are calculable in perturbative QCD, the 
^-dependence of the PDF is not and must be derived from experimental data. 

Handling the collinear and infrared divergences of the F(0(a\)) requires careful application of the factorization 
formula as illustrated in the ACOT paper. The complete 1 F contribution in VFS is given by 

l F v = F v (0(aD) = lu F v - ls F v = f N ® F^(j*g - Q t Q t ) - f N ® /f « ® F W (^Qi - ft) 

(5) 

lu F v is the 'unsubtracted' 0{a\) contribution, meaning the collinear singularities are still present, and ls F v is the 
'subtraction' term which cancels the singularities and renders 1 F finite. The V subscript indicates the use of the 
VFS. The reason for separating 1 F into these two pieces will become clear in Section 5, where the K- factors are 
discussed. Here, ^ is the perturbative 0{a\) splitting function for the process g — > Qi and is given by 

= ^r c « ln (^) w-)^ 2 -m|), (6) 

where P g ^ q (x) = t F (1 - 2.x + 2a; 2 ). In SU(N), C R = (N 2 - 1)/2N and t F = 1/2. Note that / 3 Q ' (1) = when 
fj, < Mq, and evolves continuously from zero for /x > Mq. 

4 LRSN Calculation 

The LRSN calculation uses a fixed flavor scheme (FFS), which sets the number of active flavors to a constant, 
regardless of the energy scale /z p h y of the production process. For example, when considering b production the 
number of light flavors would be four. 

The lowest-order FFS 0{a\) term is given by 

1 F F . = F P {0[a\)) - /^®F«( 7 * 5 ^Q i Q i ). (7) 

Mass factorization is not necessary at 0(a\) in FFS since the mass of the heavy quark is explicitly assumed to be 
non-zero throughout the calculation, and the initial-state heavy quark contribution is not included. Note the 1 F F 
contribution in this scheme corresponds to the 'unsubtracted' lu F v contribution of VFS. We make use of this fact 
when we compare the K- factors in Section 8. 
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LRSN has computed the Feynman diagrams shown in Figure^ to give the complete 0(a 2 s ) result. After cancelling 
the infrared divergences and performing the necessary renormalization, the mass factorization is done according to 



2 F = F(0(a 2 )) = 2 "F - 2s F = f N ® F^(j*g -> Q t Q t g) + £ /*• ® F< 2 ) ( 7 *g a - Q^) - 2s F , 

a 

(8) 

where all the contributions from initial-state partons light relative to the produced heavy quark are summed. As 
before, we recognize the first term as the 'unsubtracted' 2u F(0(a 2 )) contribution, and the second term as the 
'subtraction' term 2s F(0(a 2 )) which removes the collinear mass singularity from the radiation of a gluon by the 
incoming massless quark or gluon. Note the 'subtraction' term does not subtract the Mq — > singularity. Mass 
factorization is not necessary as Mq/Q 2 is explicitly assumed to be non-zero and the initial-state heavy quark 
contribution is assumed absent. For additional details, see references [|| and (t|. At 0(a 2 ) the FFS subscript is 
superfluous since no confusion is present at this order. 



5 K-factors 

The two calculations each include a different subset of the complete set of higher-order corrections. The comparison 
of the two approaches will determine to what extent these calculations pick up similar higher-order contributions, 
and can be used to estimate the magnitude of the corrections not included by either approach. To facilitate this 
comparison, we shall focus on the F2 K-factors of each calculation as a function of x and Q 2 . For the LRSN 
calculation, we define the K-factor to be 

1 F F + 2 F _ lu F v + 2 F 
K LR sn - — ^ = — , (9J 

r F r V 

where the Klrsn indicates the K-factor is coming from the LRSN 0(a 2 ) contributions. 
For ACOT, we define the K-factor as 



°F V + 1 F v _ lu F v + (°Fy - ls F v ) 
Kacot = — j^; = , (10) 



where Kacot indicates the K-factor is coming from the ACOT 0{oP s ) contribution. Kacot is defined in this manner 
because it is the 1 F F = lu F v contribution which is common to both calculations. Therefore, we use this common 
term to set the scale of comparison in the denominator of the K-factor. Secondly, in the threshold region (Q ~ Mq), 
the production cross section is dominated by the 0(a).) process. Viewing the 0{or s ) process as a perturbation on 
the dominant 0(a).) process, we see Equation ( fTo| ) is the natural definition of the K-factor. 

Both of these points become clear when the second form of Equation ( M ) is compared with Equation (||) . The 
higher-order contributions for the LRSN calculation are contained in the 2 F term, whereas °F — ls F yields the 
contribution for the ACOT calculation. 
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6 Parton Distributions for Heavy Quarks: The VFS Approach 

The VFS calculation of ACOT uses the CWZ renormalization scheme to incorporate the heavy quark into the QPM. 
The fundamental idea here is that at low energy scales (Q <C Mq), we do not want to treat the heavy quark as a 
constituent of the hadron. However, at high energy scales (Q ^> Mq), the mass of the 'heavy' quark is negligible 
and we should therefore treat this quark on the same footing as the other massless partons. The concept that there 
should be a democracy among the quarks in the limit Q/Mq — ■* oo forces us to introduce the heavy quark as a parton 
at some intermediate energy scale which is typically related to the quark mass. The choice of this scale is intimately 
connected to the choice of renormalization scheme Q , || . 

The motivation to include the heavy quark as a constituent of the hadron is more than aesthetic. When computing 
physical cross sections in the context of perturbation theory, we find that our perturbation expansion is not simply 
an expansion in powers of a s . As we compute higher-order subprocesses, we gain logarithms involving the various 
energy scales in the problem (such as hx(Q /Mq)). Therefore, we find that our perturbation series is actually an 
expansion in a s \n(Q / Mq). In the limit Q/Mq — ► oo, we will clearly have a divergent series unless we can resum 
these logarithmic terms. 

The QCD evolution of the parton distribution functions does precisely this resummation for the partons. It sums 
an infinite set of quark and gluon splittings, thereby taking all such logarithmic contributions into account. To be 
perfectly clear, these resummations are done in both LRSN and ACOT for the light quarks and the gluon. The only 
difference is in the way the heavy quark is treated. For ACOT, it is incorporated into the PDF, and for LRSN, it 
only enters in the final state as a product of an interaction of a light quark or gluon with the virtual photon. 

The goal of the ACOT calculation is to use this same technique to resum the multiple splittings arising from the 
QCD evolution of the heavy-quark distribution. The result is to reorganize the perturbation expansion such that the 
singular a s In (Q/Mq) terms are split into an a s \n(Q/ fj,) term which remains in the hard scattering Wilson coefficient 
and an a s Ih(ji/Mq) term which is absorbed into the heavy-quark PDF. The end result is that the perturbation series 
then becomes an expansion in powers of a s ln(Q//i), and we retain the freedom to adjust fi to optimize the convergence 
of the series. 

This is clearly the correct approach in the asymptotic limit Q/Mq — > oo but how well does it work in the 
kinematic regime of present accelerators? The only way to answer this question is to compare ACOT with a separate 
calculation such as LRSN. 

7 Kinematic Regimes 

7.1 Collinear/On-Shell vs. Large P T / Virtual 

The VFS approach (ACOT) resums an infinite subset of diagrams (Figure [|) within the PDF of the heavy quark, c/., 
via the flavor excitation (FE) (Figure ||a). Because the heavy quark is treated as a parton for this set of contributions, 
the phase space of the heavy quark is restricted to the on-shell collinear region. In the asymptotic limit \i — > oo, the 
on-shell collinear region is the dominant region of phase space but how large is this contribution in the threshold 
region fj, ~ Mq? 
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Figure 5: An illustrative subset of the class of diagrams which are summed by the QCD evolution of the heavy 
quark PDF in the variable flavor scheme (VFS) calculation of ACOT. The ® indicates a convolution between the 
heavy quark PDF (lower portion of the diagram) and the 0{oP s ) process (upper portion of the diagram). 

The FFS calculation (LRSN) does not treat the heavy quark as a parton. Instead, the heavy-quark contribution 
is included explicitly into the hard scattering via the flavor-creation (FC) process (Figure ||b). In this case, the 
heavy-quark contributions are not restricted in phase space. The splitting of the initial state gluon is only included 
up to 0(a 2 s ) but covers the collinear on-shell region as well as the region where the i-channel heavy quark is far 
off-shell. We therefore expect LRSN to provide the best results when the i-channel heavy quark is highly virtual, 
but how much of the on-shell collinear region does LRSN include? 

We want to investigate how important the roles that the flavor creation (LRSN calculation) and the flavor 
excitation (ACOT calculation) processes play when in the threshold region (Q ~ Mq), in the intermediate transition 
region (Q > Mq), and in the asymptotic region (Q ^> Mq). 

7.2 Threshold Region 

In the threshold region, we do not expect the QCD evolution of the heavy quark PDF to make a significant contribu- 
tion since the heavy-quark evolution begins at scale fi ~ Mq, and ends at scale fi *~ Q not significantly larger than 
Mq. More specifically, in the threshold region, the heavy quark will typically be produced far off-shell such that the 
dominant region of phase space is the virtual region. Since the ACOT calculation picks up primarily the on-shell 
region of phase space while the LRSN calculation picks up the entire phase space, we expect the LRSN calculation 
should make the better prediction in this region. 

In fact, if the physical threshold for heavy-quark production is below the threshold for introducing the heavy 
quark into the PDF {fi < Mq), then there will be no contribution from the heavy-quark QCD evolution since 
/q(/i < Mq) = 0. In the limit, both the °F and ls F contributions will vanish so that only lu F remains. It is 
now evident why we defined the K-factor for the ACOT calculation by Equation ([l(]). In this limit, the K-factor 
is unity, and there is no contribution from higher orders in this kinematic regimej^j On the other hand, since the 

" Note that by a choice of factorization scale, we can always ensure (although sometimes artificially) that the threshold for heavy 
quarks in the PDF's is always less than the physical heavy-quark production threshold. We will discuss the factorization scheme and the 
scale dependence in later sections. 
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LRSN calculation computes the higher order corrections for a virtual heavy-quark exchange, the result is a non-trivial 
K-factor near threshold. 

We conclude that near threshold, the dominant region of phase space is the virtual region. We therefore expect 
the LRSN calculation to determine more accurately the higher-order contributions in this region. 

7.3 Asymptotic Region 

We now consider the asymptotic region where fJ- p hy 3> Mq. In this region, we can essentially neglect the mass of 
the heavy quark in comparison to the characteristic energy scale of the process, /x p hy- Because the mass of the 
'heavy' quark is small relative to /iphy, the 'heavy' quark is easily produced on-shcll with relatively little transverse 
momentum (of order Mq) as compared to its longitudinal momentum (of order /i p hy) in the 7*-hadron center of 
momentum frame. 

Equivalently, the QCD evolution will now be significant because there is a large region (from Mq to /i) over 
which the evolution can build up the heavy-quark parton distribution. Therefore, we expect the dominant region of 
phase space is the collinear region, and that the ACOT approach should correctly resum the important higher order 
contributions. 

In this asymptotic limit, the collinear portion of 2 F in LRSN (see Figure [|) will be contained within the ACOT 
result as are the other higher-order gluon ladder graphs. If the phase space is dominated by the collinear region then 
the ACOT calculation is well suited to make accurate predictions in the asymptotic region since the dominant terms, 
the infinite set of recursive quark and gluon splittings, have been resummed. The singular a s \ti(Q/Mq) present in 
the LRSN calculation have been reorganized in ACOT (absorbed into the heavy-quark PDF) to leave only terms of 
order a s \n(Q/fi). 

If it were to happen that the terms of order 0(a^) and higher were a negligible contribution to the 0(a\) and 
0(a 2 ) terms computed in the LRSN calculation, then we would expect the ACOT and LRSN calculations to match 
in this asymptotic region. Conversely, the difference between these calculations in this region is indicative of the size 
of the terms of order 0(a^) and higher which have been resummed in the QCD evolution of the heavy quark parton 
distributions. 

While the ACOT calculation is expected to provide the more accurate results at asymptotic Q 2 , the question 
arises: what Q 2 qualifies as asymptotic? As Q 2 gets large then a s \n(Q 2 /Mq) can grow to 0(1) and spoil the 
convergence of the perturbation series. Investigating this problem in a cursory way, we present the expression for a s 
at the two-loop level to be 



where (3\ and 02 are given by 



4tt 



/?i In /i 2 / A 2 



QCD 



1 



/3 2 In In (m7A 2 



QCD) 



/??ln M 2 /A 2 



QCD 



(11) 



Pi = Y CA ~t n f tF ' (12) 

34 20 

02 = —C 2 A -—C A n s t F -AC F n s t F , (13) 
n f is the number of light flavors, C A = N, C F — (N 2 ~ l)/2N, t F = 1/2 for QCD, and N is from the SU(N) nature 
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of QCD. We see from this expression for a s that while ln(Q/Mq) is growing with increasing Q, this effect will be 
offset by the diminishing of a s with the dominant ln(/i/A) term in the denominator since /i must be chosen as a 
function of Q and/or Mq. Thus it is not clear what Q we can consider to be asymptotic, and we need to analyze 
the results of this comparison to draw definite conclusions. 

7.4 Intermediate Region 

The intermediate region is the most interesting as the heavy-quark production process is a complex interplay of all 
available mechanisms. Complications arise because we have no limiting behavior to guide us as we analyze the results 
of each calculation. What we do know is that the LRSN calculation should provide the most accurate results in the 
threshold region and ACOT should in the asymptotic. What remains to be seen is how and where the transition 
occurs from one to the other. We have little intuition about this region, and must rely on the comparison of these 
two calculations to give us insight into the physics. 

One portent of future progress in heavy-flavor production, the merging of the two calculations would be most 
effective in this region, as neither the flavor excitation process nor the flavor creation process should be the dominant 
mechanism. The three-order calculation that combines the ACOT and LRSN calculations into one should have the 
virtues of both approaches, have considerably reduced mass factorization scale dependence, and allow us to make 



very accurate predictions to compare with the data from HERA. This effort is currently underway 1 1 



8 Comparison 

We present our results using the CTEQ2 PDF set O], which begins the charm quark QCD evolution fiQ — m c = 1.6 
GeV, and the bottom quark QCD evolution at /i — m b — 5.0 GeV, For /i < (mq, the heavy-quark density in the 
hadron vanishes, Jq{x, /i < Mq) = 0. For the scale /x, we make the choice 

M \M Q + cQ 2 (l-^r ifQ 2 >M| [ ] 

with c = 1/2 and n — 2 as discussed in ACOT. 

8.1 C-Quark Production 
8.1.1 Charm F 2 vs. {x,Q 2 } 

We now present a series of plots showing the F 2 's for various values of Q 2 vs. x. The F 2 K-factors are shown to 
facilitate comparison between the different curves. The values of the F 2 structure functions are also shown on a 
log-log plot to indicate difference in an absolute sense, and to gauge the contribution to the cross section which is 
proportional to the integral of F 2 over x. 

We refer to the curves as FE for the C(a°) flavor excitation process, FC for the 0{a 1 s ) flavor creation process, 
ACOT for the complete VFS calculation and LRSN for the complete FFS calculation.^] The FC curve is not shown 
as this is the denominator in the definition of the K-factors. (Trivially, Kpc = !•) 

II Although we made a thorough comparison throughout the kinematic regime, for sake of space, only characteristic plots are shown. 
The {x, Q 2 } range was investigated in uniform logarithmic steps, which explains the odd choice of Q values presented. 
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Q = 1.778 GeV 




Figure 6: |]a) the K-factor, and ||b) the structure function F?(x, Q 2 ) for c-production at Q — 1.778GeV. 



Q = 5.623 GeV 




10" 4 10" 3 10~ 2 10" 1 1 



Figure 7: 0a) the K-factor, and 0b) the structure function F<i{x, Q 2 ) for c-production at Q — 5.623GeV. 

For Q — 1.778 GeV, (Figure ||) the FE K-factor essentially vanishes as compared to the other contributions 
because f c (x,[i) ~ 0. This result matches our expectation that "heavy-quark" partons should not contribute at 
low energy scales. The ACOT K-factor is approximately equal to one; that is, the ACOT reduces to the FC result 
in the low energy limit. This is because we have f c (x,fj,) ~ fjy(x,n) <X> fg^~\x,fj,) in the threshold region causing 
(°F V — ls F v ) ~ 0, and therefore Kacot — 1- The rise in the FE and ACOT K-factors is due to the difference 
between the one-particle and two-particle phase space factors in the threshold region (x — > x max ). The one-particle 
final state requires the partonic center-of-mass energy y/s > Mq and the two-particle final state demands y/s > 2Mq. 
This causes the K-factors for ACOT to diverge for 2M Q > y/s > Mq, where s = Q 2 (l - x)/x. The rise in the LRSN 
K-factor is attributed to the dominance of the 0(a 2 ) contributions over the 0{a\) result, as can be seen at the 
partonic level in Figures 6a and 7a of fl. The effect on the absolute value of Fi and the cross sections is negligible 
as can be seen in Figure ^|b. Clearly, in this low Q-region, the LRSN calculation is the most appropriate because the 
FC photon-gluon fusion process is the dominant one, the 0(a 2 ) corrections to that channel are included in the LRSN 
calculation and provide a non-trivial K-factor. The corrections at large x are primarily a result of the Coulomb 
exchange of a gluon between the final-state heavy quark and antiquark as well as initial-state-gluon bremsstrahlung. 

At Q — 5.623 GeV (Figure 0), the FE Fi result has increased significantly and approaches the other curves. This 
very fast evolution of the FE result makes it unreliable for predicting heavy-quark production.^] The ACOT K-factor 

** Note that for larger values of Q 2 , the FC process appears to match roughly the LRSN, and to match very well with the ACOT 
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Q = 17.78 GeV 




Figure 8: ||a) the K-factor, and ||b) the structure function F?(x, Q 2 ) for c-production at Q — \7.THGeV. 

now deviates from unity as the imperfect cancellation between °F V and ls F v signals the existence of non-trivial 
contributions from the heavy-quark PDF evolution. The fast evolution of the heavy quark in the threshold region 
(due to abundant gluons) generates important contributions for relatively low values of Q. However the subtraction 
prescription ensures the result is reliable (in contrast to the FE process), as we shall confirm when we examine the 
/z-dependence. The LRSN K-factor decreases at small x, and flattens slightly. This decrease at small x is due in 
large part to the cancellation between the mass factorization logarithmic terms and the mass factorization scale 
independent terms at large parton-photon center-of-mass energies s. For additional discussions, see pp. 192-196 of 
reference ||. 

At Q = 17.78 GeV (Figure §), the FE, ACOT, and LRSN results have similar shapes. The K-factors are 
all monotonically increasing vs. x, with a large rise in the threshold (x — > x max ) region. In the range above 
Q = 17.78 GeV, the general characteristics are similar to this Figure ||. 

8.1.2 Charm F 2 vs. fj, 

We now present a series of plots of the /z-dependence of F 2 for particular x and Q 2 values accessible to HERA. 

For x = 0.001 and Q = 3.162 GeV (Figure ||), the FE result has a large /i-dependence. The FE process is driven 
by the heavy-quark PDF and vanishes at small y, because / c (x, \i) = for fi < m c . As fi increases above m c , / c (x, fi) 
increases quickly due to the abundance of gluons and the dearth of heavy quarks near threshold. Clearly, the FE 
result is extremely sensitive to the particular choice of scale, as well as the choice of factorization scheme which 
determines where to start the heavy-quark evolution. As anticipated, the apparent agreement between the FE result 
and the other results as seen in the previous subsection is merely an accident due to a prudent choice of scale, and 
cannot provide stable results. The LRSN result is quite flat, and the ACOT calculation exhibits a rather marked 
dependence on the scale choice, although substantially less than the FE result. 

For x = 0.001 and Q = 10 GeV (Figure |l0|), the ACOT and the LRSN calculations exhibit a comparable and 
small dependence on the scale choice. In this region, neither is a clear improvement over the FC result. 

For x = 0.01 and Q = 10 GeV (Figure |ll|), both the LRSN and ACOT calculations are essentially flat, and both 

calculation. This agreement is fortuitous and depends on a judicious choice for the factorization scale fi. As we shall see shortly, the FC 
result is very sensitive to jx. 
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Figure 9: The /z-dependence of F2(x,Q 2 ) with {x, Q} — {10 3 , 3.162 GeV} for c-quark production. 
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Figure 10: The /x-dependence of -F 2 (x, Q 2 ) with {x, Q} = {10 3 , 10 GeV} for c-quark production. 

appear to be an improvement over the FC result. 

For x = 0.1 and Q = 31.62 GeV (Figure |l2]) , the ACOT calculation exhibits slightly less dependence than the 
LRSN result, and both are considerably less scale-dependent than the FC and FE results. 

The general trend revealed here confirms the reduced /i-dependence of the ACOT and LRSN results over their 
one-order counterparts. Furthermore, we also have confirmation that LRSN is the appropriate calculation near 
threshold and ACOT in the asymptotic region, as shown when Q 3> Mq. 

8.2 B-Quark Production 
8.2.1 Bottom F 2 vs. {x,Q 2 } 

We have also compared the ACOT and the LRSN results for b-quark production. We find very similar features in 
Figure |l3| with Q = 5.623 GeV for b production as we did in Figure ^ with Q — 1.778 GeV for c. This is as expected 
since the physics is set by the relevant mass ratio in the process: Mq/[i. The LRSN calculation has a non-trivial 
K-factor, the FE result essentially vanishes, and the ACOT K-factor is essentially unity. 

For Q = 10 GeV (Figure |l4|), the FE results are roughly comparable to the other curves as the heavy-quark 
evolution turns on. The ACOT K-factor deviates from the FC result as the difference between °F V and ls F v again 
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Figure 11: The //-dependence of F%{x, Q 2 ) with {x, Q} = {10 2 , 10 GeV} for c-quark production. 
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Figure 12: The //-dependence of F 2 (x, Q 2 ) with {x, Q} = {10 _1 , 31.62 GeV} for c-quark production. 

signals the existence of non-trivial contributions from the heavy-quark PDF evolution. The LRSN K-factor decreases 
at small- a;, and flattens slightly again for the reasons discussed concerning Figure 0. 

In Figure [l5|, we display the results for Q = 31.62 GeV. The ACOT exhibits a significant difference from the 
FC process as the heavy-quark evolution continues. As before, the FE, ACOT, and LRSN results have similar 
shapes, and each K-factor is monotonically increasing as x increases. In the range above Q = 31.62 GeV, the general 
characteristics are similar to this figure. 

8.2.2 Bottom F 2 vs. [i 

We now examine the //-dependence of b production at various points in x and Q 2 space relevant to the region 
accessible at HERA. 

For x — 0.001 and Q = 10 GeV (Figure |l6]), the FE result again has a large /i-dependence because this result is 
closely tied to fb(x : f£). The ACOT result is a clear improvement over the FE result, but the FC result exhibits still 
less dependence and is comparable with the LRSN calculation. 

For x — 0.01 and Q = 10 GeV (Figure [l7|), the FE result still has a large //-dependence, and both the ACOT and 
the LRSN results improve upon the FC result. 
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Figure 13: |l3|a) the K-factor, and |T^b) the structure function F2(x, Q 2 ) for 6-production at Q = 5.623GeV. 
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Figure 14: |l4|a) the K-factor, and p^|b) the structure function F2(x,Q 2 ) for fo-production at Q = WGeV. 

For x — 0.1 and Q = 31.62 GeV (Figure |l8|), the LRSN and the ACOT calculations are comparable and both 
show substantially less /^-dependence than the FE or FC results. 

We can make some general observations regarding the /i-dependence for both charm and bottom production. 
For all values of Q, the FE process is increasing with fi due to the increasing heavy-quark PDF. In contrast, the 
0(al) FC process (driven by gluons) is decreasing with /t largely due to the decrease in a s (/i). The two-order 
calculations (ACOT and LRSN) that have compensating contributions to cancel out some of the ^-dependence. 
Specifically, ACOT combines pieces of the FE and FC processes (together with a subtraction term) to yield a result 
that has substantially less /x-dependence than either result in the large Q region. LRSN effectively has the 0(a 2 ) 
FE contribution as the collinear heavy-quark part of phase space is included, negating some of the /(-dependence of 
the FC channel. 

We note that in general, the ACOT result exhibits minimal //-dependence at larger values of Q and x. At 
large Q, we are closer to the asymptotic region where the VFS approach is expected to be superior. The increased 
/i-dependence at low x arises mainly from the quark-initated contributions. At a given order, we expect the quark- 
initiated contributions to be less than the gluon-initiated ones because f q (x,[i) <C f g (x,fi) B. However, because 

r<2) (2) 

the second-order Altarclli-Parisi splitting kernels PqJ, q (x) and Pg^ q {x) contain singular 1/x terms whereas the first- 
order Pq}l q (x) and Pg\}, q {x) do not, the evolution of the quark distribution at small-x will be completely dominated 
by the second-order kernel rather than the first-order kernel p2|. Consequently, contributions from higher order 
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Figure 15: |l5|a) the K-factor, and |l5|b) the structure function F2(x, Q ) for 6-production at Q — 3l.62GeV. 
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Figure 16: The /i-dependence of F 2 (x,Q 2 ) with {x, Q} = {10 3 , 10 GeV} for 6-quark production. 

quark-initiated processes should cancel the above /^-dependence. This work is in progress pc| |. 

In examining the /i-dependence plots, the reader may have noticed for some x and Q 2 values chosen, particularly 
at small x, the results of ACOT and LRSN never crossed. The conventional wisdom says the /^-dependence represents 
the uncertainty from uncalculated higher-order contributions. What are the implications when the /i-dependence 
curves do not cross when a significant region along the /t-axis has been traversed? At lower values of Q 2 , the difference 
between ACOT and LRSN is expected due to the limited evolution of the heavy quark PDF. At asymptotic values 
of Q 2 , the difference between ACOT and LRSN is within the range suggested by the /* dependence. However, in the 
intermediate range, the range of the /i-dependence may underestimate the full theoretical uncertainty. 



8.2.3 Total Charm and Bottom Cross Section 

To estimate the theoretical uncertainty for the total cross section, we have computed 



F 2 (Q 2 



io- 4 



dxF 2 (x,Q 2 ) 



(15) 



which is proportional to the dominant contribution for the heavy-quark production cross section and show the ratio 



of F 2 COT I and F^ RSN jF^ c for c production in Figure 19. For low Q, the differences are on the order of 30%. 
As Q increases, however, the differences decrease to ~ 12 ~ 15% as Q ~ 100 GeV. Similar behavior occurs for b 
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Figure 17: The ^-dependence of F 2 (x, Q 2 ) with {x, Q} — {10 2 , 10 GeV} for 6-quark production. 
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Figure 18: The ^-dependence of F 2 (x,Q 2 ) with {x,Q} = {10 1 , 31.62 GeV} for 6-quark production. 

production, as seen in Figure [|(]. Some of this difference may possibly be attributed to choice of scale or PDF set. 
However variations in scale or PDF set should not be solely responsible for this difference. We may have to await 
the data from HERA to resolve this discrepancy. 



9 Conclusions 

We have outlined the strengths and weaknesses of both the VFS (ACOT) and the FFS (LRSN) calculation. We 
summarize the highlights below. 

• In the threshold region, the FFS (LRSN) calculation yields the most stable and reliable results due to the 
domination of flavor creation. 

• In the asymptotic region, the VFS (ACOT) calculation provides the best results because of the dominance of 
the collinear heavy-quark contribution. 

• In the kinematic range spanned by HERA, the a s \og[Q 2 /Mq\ terms are not a significant problem for the FFS 
(LRSN) calculation. 
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Figure 19: The structure function ^(Q 2 ) scaled by the flavor creation (FC) ^(Q 2 ) for c-quark production. 
F%(x, Q 2 ) has been integrated over x from x — 10 -4 to x max as defined in Equation (IT 
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Figure 20: The structure function ^(Q 2 ) scaled by the flavor creation (FC) i^Q 2 ) for &-quark production. 
F2(x, Q 2 ) has been integrated over x from x — 10 -4 to x max as defined in Equation (II 



• In the VFS (ACOT) calculations, the heavy-quark PDF's can yield significant contributions at relatively small 
scales, (i.e. (j./Mq ~ 3). 

• While the flavor excitation (FE) process can closely match the two-order results with a judicious choice of the 
scale fi, the large scale dependence makes this unreliable in computing structure functions. 

• Likewise, while the flavor creation (FC) process is a good starting point in the threshold region. However, the 
LRSN calculation indicates that the corrections to this naive estimate can be as large as 100%. 

The range of x and Q 2 we have presented reflect the region accessible at HERA. We note that the difference 
between the LRSN and ACOT calculations above threshold is suggestive of higher order contributions yet to be 
included. As such, the results of this comparison indicate that a combining of the LRSN and ACOT calculations 
in a consistent fashion (with the additional mass factorizations required) should allow us to make predictions based 
upon a three-order result that combines the best attributes of both calculations. The result should be a calculation 
that will provide an important test of perturbative QCD when compared with the results from HERA. 



18 



Acknowledgements 

The authors would like to thank John Collins, Jack Smith, Davison Soper and Wu-Ki Tung for useful discussions. 
This work is partially supported by the U.S. Department of Energy Contract No. DE-FG05-92ER-40722, by the 
Texas National Research Laboratory Commission, and by the Lightner-Sams Foundation. F.O. is supported in part 
by an SSC Fellowship. 

10 Bibliography 
References 

[1] A.C. Bcnvcnuti et al, (BCDMS Collaboration), Phys. Letts. B223 (1989) 485; B237 (1990) 592; P. Berge et 
al, (CDHSW Collaboration), Z. Phys. C35 (1987) 443; P.Z. Quintas et al, CCFR collaboration, Phys. Rev. 
Letts. 71 (1993) 1307; W.C. Leung, et al, (CCFR Collaboration), Phys. Letts. B317 (1993) 655; J.J. Aubert 
et al, (EMC Collaboration), Nucl. Phys. B259 (1985) 189; B272 (1986) 158; P. Amaudruz et al, (NMC 
Collaboration), Phys. Letts. B295 (1992) 159; A. Bodek et al, (SLAC-MIT Collaboration)! 3 /^. Rev. Letts. 50 
(1983) 1431; 51 (1983) 534. 

[2] A. Bazarko, et al, (CCFR Collaboration), Nevis R#1502. 

[3] G.A. Schuler, F.I. Olness, J. Blumlein, and W.-K. Tung in Proceedings of the 1990 DPF Summer Study on High 
Energy Physics: Research Directions for the Decade, Snowmass, CO, p. 152, (1992). 

[4] M.A.G. Aivazis, F.I. Olness, and W.-K. Tung, SMU-HEP/93-16, to appear in Phys. Rev. D50, (1994); 

M.A.G. Aivazis, J.C. Collins, F.I. Olness and W.-K. Tung, SMU-HEP/93-17, to appear in Phys. Rev. D50, 
(1994). 

[5] J.C. Collins, F. Wilczek, and A. Zee, Phys. Rev. D18, 242 (1978). 

[6] E. Laencn, S. Riemersma, J. Smith, and W.L. van Necrvcn, Nucl. Phys. B392, 162 (1993). 
[7] E. Laencn, S. Riemersma, J. Smith, and W.L. van Necrvcn, Nucl. Phys. B392, 229 (1993). 
[8] J.C. Collins, W.-K. Tung Nucl. Phys. B278, 934 (1986). 

[9] Sijin Qian, The CWZ subtraction scheme (A new renormalization prescription for QCD) and its applications, 
Illinois Inst, of Tech. Ph.D. thesis, UMI-85-17585-mc (microfiche), May 1985. 

[10] P. Agrawal, F.I. Olness, S. Riemersma, and W.-K. Tung, SMU preprint in preparation. 

[11] J. Botts, H.L. Lai, J.C. Morfm, J.F. Owens, J. Qiu, W.-K. Tung and H. Weerts, Version 2 CTEQ Distribution 
Functions in Parametrized From, MSU preprint in preparation. 

[12] W.-K. Tung, Nucl. Phys. B315, 378 (1989). 



19 



